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We propose a novel mult i- component system of nonlinear equations that generalizes 
the short pulse (SP) equation describing the propagation of ultra-short pulses in optical 
fibers. By means of the bilinear formalism combined with a hodograph transformation, 
we obtain its multi-soliton solutions in the form of a parametric representation. Notably, 
unlike the determinantal solutions of the SP equation, the proposed system is found to 
exhibit solutions expressed in terms of pfaffians. The proof of the solutions is performed 
within the framework of an elementary theory of determinants. The reduced 2-component 
system deserves a special consideration. In particular, we show by establishing a Lax pair 
that the system is completely integrable. The properties of solutions such as loop solitons 
and breathers are investigated in detail, confirming their solitonic behavior. A variant of 
the 2-component system is also discussed with its multisoliton solutions. 
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I. INTRODUCTION 

The short pulse (SP) equation was derived as a model nonlinear equation describing the 
propagation of ultra-short pulses in isotropic optical fibers.^ We write it in an appropriate 
dimensionless form as 

Uxt = U + -{U^)^^, (1.1) 

where u — u{x, t) represents the magnitude of the electric field and subscripts x and t 
appended to u denote partial differentiations. The SP equation has appeared for the 
first time in an attempt to construct integrable differential equations associated with 
pseudospherical surfaces.^ The integrability, soliton solutions and other features of the SP 
equation common to the completely integrable partial differential equations (PDEs) have 
been studied from various points of view.^~^° See also Ref. 11 for a recent review article 
on the SP equation which is mainly concerned with soliton and periodic solutions and 
their properties. It also provides a novel method for constructing multiperiodic solutions 
by means of the bilinear transformation method. 

There exist a few generalizations of the SP equation to the two-component systems that 
take into account the effects of polarization and nonisotropy. One is due to Pietrzyk et al. 
They proposed the following three integrable vector (or two-component) SP equations:^^ 



Uxt^U+ -^{U^ + ^Uv'^)xx, Vxt^V + -^{v^ + ^u'^v)xx, (1-2) 

U^t^U+ ^{U^ - 3uv'^)a:x, Vxt ^ V - ^{v^ " 3^^^)^;^;, (1.3) 

Mxt = « + ^{u^)xx, v^t = v + ^{u'^v)^^. (1.4) 
Another one is given by Sakovich:^^ 

1 1 

Uxt^u+ -{u^ + uv'^)^^, v^t ^ V + -{v^ + u'^v)^^, (1.5) 

u^t^u + ^{u^)xx, v^t^v + hu'^v)^^. (1.6) 



As pointed out by Sakovich/^ the two systems (1.2) and (1.3) can be reduced to the SP 
equation (1.1) by appropriate dependent variable transformations. Indeed, introducing 
the new variables p and q hj p = u + v,q = u — v, the system of equations (1.2) can 
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be decoupled and both p and q satisfy the SP equation (1.1) while for (1.3), the trans- 
formation p — u -\- \v and q — u — iv leads to the two decoupled SP equations as well. 
On the other hand, the integrability of the systems (1.5) and (1.6) was investigated by 
means of the Painleve analysis. Sakovich showed that the above two systems pass the 
Painleve test, providing a strong indication of their integrability. Nevertheless, their Lax 
representations, conservations laws and soliton solutions have not been obtained as yet 
for the systems. 

The purpose of this paper is to propose a novel multi-component analog of the SP 
equation and construct its multisoliton solutions. The system of equations presented here 
is composed of the following coupled nonlinear PDEs for the n variables Ui{i — 1,2, ...,n): 

Ui,xt^Ui + ]^{Fui^^)^, (i = 1,2, ...,n), (1.7a) 

with 

-F= ^ CjkUjUk. (1.76) 
i<i<fc<" 

Here, Cjk are arbitrary constants with the symmetry Cjk = Ckj{j,k = 1,2, ...,n). For the 
special case of n = 2 with Ci2 = 1, this system becomes 

1 1 
Uxt^u+ -{uvux)x, Vxt = v + -{uvvx)x, (1.8) 

where u — ui and v — U2- Obviously, if we put u — v, then (1.8) reduces to the SP 
equation (1.1). A simple transformation recasts (1.8) to the system of equations 

li V = 0, then this system reduces to the SP equation (1.1). The present paper is 
organized as follows. In Sec. II, we summarize an exact method of solution for the 
SP equation which will be suitable for application to the multi-component system. In 
Sec. Ill, We show by applying the standard procedure of the bilinear method that the 
system of equations (1.7) can be transformed to a coupled system of bilinear equations and 
obtain the multisoliton soliton solution in the parametric form. Notably, the tau-functions 
constituting the solution are expressed in terms of pfaffians unlike the determinantal 
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solutions of the SP equation.^ The proof of the multisohton solution is, however, performed 
with use of an elementary theory of determinants without recourse to the pfaffian theory. 
In Section IV, we consider the system (1.8). In particular, we demonstrate that it is a 
completely integrable system by establishing a Lax pair. The multisoliton solution to 

the system is reduced from that of the n-component system. The properties of the 1- 
and 2-soliton solutions will be investigated in detail. Subsequently, wc briefly discuss the 
system (1.9). In Sec V, we conclude this study with a short summary and discuss some 
open problems associated with the multi-component SP equations. 

II. SUMMARY OF THE EXACT METHOD OF SOLUTION 

Here, we give a short summary of the exact method of solution for the SP equation. 
Although we have employed some nonlinear transformations to reduce the SP to the 
integrable sine-Gordon (sG) equations,^'^'^° we provide a different approach which is more 
suitable for solving the system of equations (1.7). 

A. Hodograph transformation 

We first introduce the hodograph transformation {x, t) {y, r) by 

dy = rdx -\- -u^rdt, dr = dt, (2.1a) 

where r(> 0) is a function of u to be determined later. Using (2.1a), the x and t derivatives 

are rewritten as 

d d d d 1 o d /^^,\ 

It follows from (2.1b) that x — x{y, r) satisfies the system of hnear PDEs 

1 

~ ~' — (2-2) 
Equation (1.1) is then transformed into the form 

Uyr = XyU. (2.3) 



The form of r can be determined by the solvability condition of the system (2.2), i.e.. 



XyT — XT-y. Indeed, this immediately gives = uUyr"^. On the other hand, it follows 
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from (2.2) and (2.3) that u — rUyr- Eliminating the variable u from both relations, one 
has rr — UyUyrT^. If we impose the boundary conditions u{±oo,t) — 0,r(±oo,T) = 1, 
then we obtain r'^ — (1 — Uy)~^ after integrating this relation with respect to r. Since 
Uy = Ux/r by (2.1b), we can rewrite this expression into the form 



(2.4) 



The above relation has been used to transform the SP equation into the form of conserva- 
tion law Tt — (w^r/2)^. If one introduces a new variable by = sin 0, then satisfies 
the sG equation (f)yT — sin 0. This equation was the starting point in constructing mul- 
tisohton solutions of the SP equation.^ Below, we develop an alternative method using 
(2.3) which will be relevant to application to the multi-component system. 

B. Pcirametric representation of soliton solutions 

The sohton solutions of Eq. (2.3) are constructed by a direct method using the bilinear 
formalism. To this end, we first introduce the following dependent variable transforma- 
tions for u and x 

9 



u 



r 



h 



(2.5) 



(2.6) 



where /, g and h are tau- functions. Note that we may add an arbitrary constant on the 
right-hand side of (2.6), if necessary. The second equation of (2.2) is then transformed to 
the the bilinear equation 

2Drh-f + g^^0, (2.7) 
where the bihnear operators Dr and Dy are defined by 

Kn;f-,= {^-^) (^-^) nr,MrW) 



T'=T,y'=y 



(m,n = 0,l,2,...). 

(2.8) 



On the other hand, Eq. (2.3) becomes 



iJy 
P 



(2/r + h)- -^{UOy + fyOr + fyrO + ghy) + -.{Oyr - g) = 0- 



f 



f 



(2.9) 
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We can decouple Eq. (2.9) to a set of equations 



2/, + /i = 0, (2.10) 

frgy + fygr + fyrQ + Qhy - f {Oyr - Q) = 0. (2.11) 

Substituting h from (2.10) into Eqs. (2.11) and (2.7), we obtain the following system of 
bilinear equations for / and g: 

DyDJ-g^fg, (2.12) 
Dlf-f^lg'. (2.13) 
It then follows from (2.6) and (2.10) that 

(2.14) 

Thus, the soliton solutions of the SP equation are given by the parametric representation 
(2.5) and (2.14) in terms of the tau-functions / and g. In the simplest case of the 1-soliton 
solution, the solutions to Eqs. (2.12) and (2.13) are easily found to be as 

/ = l + e2«, g=-e^, ^ = py + -r + ^o, (2.15) 
P P 

where p and ai'e constants related to the amplitude and phase of the soliton, respectively. 
The corresponding parametric representation of the 1-soliton solution is derived from 
(2.5), (2.14) and (2.15). It reads 

2 2 

ii=-sech^, x = y temh. ^ + xq, (2-16) 

P P 

where Xq — — p/2. For real p and ^o, the solution takes the form of a loop soliton. 
C. Remcirk 

We have already shown that the SP equation can be transformed into the sG equation 
and obtained the parametric representation of the A^-soliton solution. Actually, it reads^ 

ii = 2i(ln£j , x = y-2(ln/'/)r, (2.17) 
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where / and /' are tau-functions for the sG equation 0^,- = sin 0, (f) — 2i ln(/'//) and 
they satisfy the bihnear equations 



The exphct forms of the tau-functions are given by 



(2.18) 



M=0,1 

f'=Yl 

/i=0,l 



' ' l<j<k<N 



i=i 

N 



.3=1 l<j<k<N 



where 



1 



= PjV + -r + ^jo, U = l,2,...,iV), 
Pj 



Pj - Pk 



{j,k = 1,2,. ..,N; J j^k). 



(2.19a) 



(2.196) 



(2.20a) 



(2.206) 



.Pj +Pk, 

Here, pj and are arbitrary complex-valued parameters satisfying the conditions pj ^ 
±Pk for j ^ k and N is an arbitrary positive integer. The notation ^^=q i implies the 
summation over all possible combination of /xi = 0, l,/i2 = 0, 1, ...jfiN = 0, 1. Thus, we 

have two different expressions for the parametric soliton solutions of the SP equation, i.e., 
one is (2.5) with (2.14) and the other is (2.17). Wc can show that the tau-functions / 
and g are related to the tau-functions / and /' by the relations 



f = f'f, g = 2iDj'-f, 
which will be inferred by comparing (2.5) and (2.14) with (2.17). 



(2.21) 



III. MULTI-COMPONENT SYSTEM 

Let us now consider the multi-component system (1.7). The procedure for obtaining 
the parametric representation of soliton solutions parallels that developed in Sec. II for the 
SP equation. Hence, we omit the detail of the derivation and write down the final results. 
Specifically, we give the parametric representation of soliton solutions and associated 
system of bilinear equations corresponding to Eqs. (2.12) and (2.13). Then, we present 
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the explicit form of the multisohton solution of the bilinear equations. Last, the proof of 
the multisoliton solution is performed by using an elementary theory of determinants. 
A. Pcirametric representation of soliton solutions 

If we use the hodograph transformation (2.1a) with F given by (1.7b) in place of 

dy — rdx + -Frdt, dr — dt, (3.1) 

we then obtain the equations corresponding to Eqs. (2.2) and (2.3) which are given 
respectively by 

_ 1 _ F 

— ~i — (3.2) 

Ui,yT = XyUi (i = 1,2, ...,n). (3.3) 

The solvabihty condition for Eqs. (3.2) gives the exphcit form of in terms of Ui^y {i — 
1, 2, n) as 

r' = 73^. ■ (3.4a) 

If we use the relation Uj^y — Uj^xl"^-! then we can rewrite (3.4a) in terms of the original 
variable Ui^^ (i = 1, 2, n) 

= 1 + CjkUj,xUk,x- (3.46) 

l<j<k<n 

The parametric representation of the soliton solutions takes the form 

Mi=y, (« = 1,2, ...,n), a; = y-2y, (3.5) 
where the tau-functions / and gi{i = 1, 2, n) satisfy the system of bilinear equations 

DyDrf ■ Qi^ fgi, (i = l,2,...,n), (3.6) 
Dlf • / = ^ Yl ^ik9j9k- (3.7) 

l<j<fc<n 

It follows from (3.2)-(3.4a) that Ui — Ui{y,T) obey a closed system of PDEs 

M„ (^ = l,2,...,n). (3.8) 



-^/l Y2l<j<k<n ^jkUj^yUk^y 
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Furthermore, if we introduce the new variables Vi by Vi — Ui^y{i — 1,2, then the 

above system can be recast to 



d_ 
dy 



B. Multisoliton solution of bilinear equations 



1,2,. 



(3.9) 



We first introduce vectors and matrices. Subsequently, we present the explicit multi- 
soliton solution of the bilinear equations (3.6) and (3.7). 
1. Definition 

Let a, b, c and be row vectors having M components 

a= (ai,a2,...,aM), b = (6i, 62, J'm), c = (ci, C2, cm), 

d^idi,d2,...,dM), = (0,0,...,0), (3.10a) 

and e^ {i = 1, 2, n) be M-component row vectors defined below 

ei = (M^,OA^,---,ei = ( ), 

Ml M-Mi Mi+-+Mi_i Mi M-{Mi+---+Mi) 

...,e„ = ( 0^0^ , M^), (3.106) 

Mi+---+M„-i M„ 

where M and Mj(i = 1,2, ...,n) are positive integers satisfying the condition Y17=i-^'i ~ 
M. 

The following types of matrices appear in the process of proving the multisoliton 
solution: 

A, 



D{aL, b; c, d) 



/ Am 


Im\ 


\-Im 


Bm) 


Im 




Bm 0^ 


0^ 











0^ 



L'(a, ei;b,ej) 



Im 










Bm 






(3.11a) 













(Am 


Im 


b^ 






—Im 


Bm 


0^ 






a 















ej 





0^ 










(3.116) 
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D{a,h,ei]C,d,ej) = 



/Am Im 

-Im Bm 0^ 0^ ej 

a 

b 

\ y 



(3.11c) 



where Am — {(iij)i<i,j<M and Bm — {bij)i<i,j<M are M x M skew-symmetic matrices, Im 
is the M X M unit matrix and the symbol T denotes the transpose. 

Pi - V-i ^ 



The element aij of the matrix Am is given by 

= -^^e^'^^^- = -^-^z,z,, = 1,2,...,M), (3.12) 

Pi + Pj Pi + Pj 

.20a) and Zi = cP. To specify the matrix Bm, let Si {i = 1, 2, n) 



Qij — 



where is defined by (2 

be n disjoint sets consisting of positive integers 



Then 



Si = {1, Ml}, S, = {M, + M2 + --- + M,_i + 1, Ml + ■ ■ ■ + M,}, 
= {Ml + M2 + • • • + M„_i + 1, Ml + • • • + M„}. 



(3.13) 



bjj,!/ — . Cij 2 2 ■ ^ ^ 



,5^, 1/ e -Sj (//, 1/ = 1, 2, M {pL^uy,iJ ^ 1, 2, n (i 7^ j)), 

(3.14) 

6^,^ = if /X and u belong to the same set and — for all Thus, Bm has the 

striiptiirp 



Bm — 



^ OmixMi BmixM2 ■■■ BmixM„^ 
~Bm^xM2 OM2XM2 ■■■ Bm2xM„ 



\ Bj^^y^Mn Bm^^Mu ■■■ OMnXMnJ 

BMixMj = {bij.i.)ij.eSi,ueSj (1 < « < J < n), 
OMixMi : MiX Mi null matrix {i = 1, 2, n). 



(3.15a) 



(3.156) 



2. Multisoliton solution 

Now, we state our main result. 

Theorem 3.1: The multisoliton solution of the system of bilinear equations (3.6) and 
( 3. 7) is given by the following form: 

f = VF, F = \D\, (3.16a) 
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Qi^wGi, = |i:)(-z, -ei;z,ej)|, (i = 1,2, (3.166) 

where z is the M-component vector z = (e^% e^^, e^^'^). The parametric solution Ui (3.5) 
constructed from these tau-functions contains Mi solitons for each i. 

Note that / and gi are pfaffians since each one of them is represented by the square 
root of the skew-symmetric determinant of even oder. This fact is in striking contrast to 
the tau-functions of the A^-sohton solution for the SP equation which can be represented 
by determinants. 

C. PROOF OF MULTISOLITON SOLUTION 

1. Bcisic formulas for determinants 

Let A — {aij)i<ij<M be an M x M matrix and Aij be the cofactor of the element 
Qij. Then, the following three formulas for determinants are employed frequently in our 
analysis 

M 

A\z - Aijaihj, (3.18) 

|^(a,b;c,d)||A| = |A(a;c)||A(b;d)| - |^(a; d)P(b; c)|. (3.19) 

The formula (3.17) is the differential rule of the determinant and (3.18) is the expansion 
formula for a bordered determinant with respect to the last row and column. The formula 
(3.19) is Jacobi's identity and it will play a central role in the proof of the multisoliton 
solution. 

2. Differential formulas 

We give various differential formulas for the determinants F and introduced in 
(3.16) which are necessary for the proof of the solution. The following formulas are 
derived easily with use of (3.17) and (3.18) as well as the relation |D(— z;z)| = which 
follows from the fact that the skew-symmetric determinant of odd order is identically zero. 
Hence, we quote only the results: 

F, = -2|L>(-z;z,)|, (3.20a) 



A a^ 
b z 
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= -2|D(-z^;z)|, (3.206) 

Fyr = -2\D{-Zr;Zy)\-2\D{-z,-Zr;z,Zy)\, (3.20c) 

Frr = -2\D{-Zrr;z)\ - 2|D(-z, -z^; z, z^) | , (3.20d) 

Gi^y = 2\D{-z,-ei;Zy,ei)\, (3.21a) 

Gi,r^2\D{-Zr,-ei;z,ei)\, (3.216) 
Gi^yr = 2|D(-z, -ei;z,ei)|+2|D(-z^, -e^; Zj^, ei)|+2|D(-z, -z^, -e^; z, z^^, ej)|, (3.21c) 
where the M-component vectors Zj^,z^ and z^-^ are given respectively by 

, . . . , /e^i e^2 e^^\ f e^^ e^^ ^^mX 

\Pi P2 pmJ \pi pi PiiJ 

(3.22) 

3. Proof of Eq. (3.6) 

First, we show that the tau-functions (3.16) for the multisohton solution satisfy the 
bihnear equation (3.6). To this end, we substitute / and Qi from (3.16) into Eq. (3.6) to 
obtain 

^ {pFyr - -j^FyF^ {GiGi,yT - -^Gi,yGi,T^ - -^{FyGi,T + FrGi^y) = FGi. (3.23) 

We compute three terms on the left-hand side of (3.23) separately. Using (3.20a)-(3.20c) 
and the relation 

|D(-z,-z,;z,z,)||D| = -|D(-z;z,)||D(-z,;z)|, (3.24) 

which follows from Jacobi's identity and the identity \D{—z; z)| = 0, the first term on the 
left-hand side of (3.23) reduces to 



G^ 
2F 



(^FFy, - \FyF^ = -\D{-z,-Zy)\Gi. (3.25) 
Next, it follows from (3.21a)-(3.21c) that 

GiGi^yr - ^Gi^yGi^r = 2\D{-z, -ef, z, e^) | ^\D{-z, -e^; z, e^) | + \D{-Zr, -e^; Zy, e,) | 
+ \D{-z, -Zr, -ei;z,Zy,ei)\^ - 2\D{-z, -ei;Zy, ei)\\D {-Zr, -ei;z,ei)\. (3.26) 
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Referring again to Jacobi's identity and the identity \D{—ei;ei)\ — 0, one has 

|D(-z, -z^, -ei;z,Zy,ei)\\D{-ei;ei)\ 

= \D(-z, -ei;z,ei)\\D(-Zr, -ei;Zy,ei)\ - \D(-Zr, -ei;z,ei)\\D(-z, -ei;Zy,ei)\ = 0, 

(3.27) 

which, introduced into (3.26), simphfies the second term on the left-hand side of (3.23) 

^ (^GiGi^yr -^Gi^yGi^r^ = f||D(-z, -z^,-ei;z,Zj/,ei)| (3.28) 

Last, the formulas (3.20a), (3.20b), (3.21a) and (3.21b) give simply the third term on 
the left-hand side of (3.23) 

-^{FyGi^r + FrGi,y) = \D{-z; Zy)\\D{-Zr, -ef, z, ei)\ + \D{-Zr;z)\\D{-z, -ei;zy,ei)\. 

(3.29) 

Substituting (3.25), (3.28) and (3.29) into (3.23), the equation to be proved becomes 
\D\\D{-z,-Zr,-ei;z,Zy,ei)\ - \D(-Zr;Zy)\\D(-z, -ei;z,ei)\ 

+ \D{-z;Zy)\\D{-Zr,-ei;z,ei)\ + |D(-z^; z)||D(-z, -6^; Zj^, e^)] = 0. (3.30) 

The following formula can be verified by applying Jacobi's identity twice to the right-hand 
side of (3.31): 



|L>(a;a')| 


l^(a;b')| 


\D{a;c')\ 






|D(b;a')| 


|D(b;b')| 


|D(b;c')| 


= |L>|2|L>(a,b,c;a',b',c')|. 


(3.31) 


|D(c;a')| 


\D{c;h')\ 


\D{c;c')\ 







Assume that \D\ ^ 0. Then, multiplying (3.30) by \D\ and using Jacobi's identity as well 
as the identities |£)(— ei;ei)| = \D{—z\z)\ — 0, the resulting relation reduces to (3.31) 
with the identification a = — z, b = — z^, c = — Cj, a' = z, b' = z^^, c' = ej. This completes 
the proof of Eq. (3.6). 
4. Proof of Eq. (3.7) 

We proceed to the proof of Eq. (3.7). By using / and gi from (3.16) and noting the 
symmetry = c^j, we transform it to the form 

1 " 

FFtt ~ -^r = ^ CjkGjGk- (3.32) 

j,k=l 
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If we substitute (3.16b), (3.20b) and (3.20d) into (3.32) and use the following relation 
with j — k 



\D{-ej;z)\\D{-ek;z)\ = |D||D(-z, -e^-; z, ejt)|, (j, A; = 1,2, ...,n), (3.33) 

which comes from Jacobi's identity, we recast (3.32) in the form 

1 " 

2|D|||D(-z;z,,)| - |D(-z,-z,;z,zO|} = ^Yl Cjk\D{-ej;z)\\D{-ek;z)\. (3.34) 

j,k=i 

Last, replacing the right-hand side of (3.34) by the right-hand side of (3.33) and dividing 

the resultant equation by 2\D\, the equation to be proved reduces to the following linear 
relation among determinants: 

1 " 

|^(-z;zrr)| - \D{-z, -z^;z,z^)| = g XI -Gj;z,ek)\. (3.35) 

j,k=i 

We now start the proof of (3.35). Define the (2M -|- 1) x (2M -|- 1) skew-symmetric 
matrix D' = (d-j)i<jj<2M+i by 

(Am Im z^\ 
-Im Bm 0^\. (3.36) 
-z / 

Let Dij and D'-j be the cof actors of the elements dij and d'^^p respectively and Dij ki and 
^ij,ki second cofactors. Expanding the cofactor D'j^_^-^f^_^^ with respect to ith row, we 
obtain 

M M 

^M+j,M+i — X^ B)j^M+j,kM+i^ik + DiM+j,kM+i ZiZk, {i, j — 1, 2, M). (3.37) 

k=l k=l 

Similarly, referring to the structure of the matrix Bm defined by (3.15), the expansions 
of Dij and D^j with respect to the (M -|- i)th column read 

M 

Dij = X DiM+kjM+ihi, = 1, 2, M), (3.38) 

k=l 
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M 

D'ij = J2^iM+k,jM+ihi, = 1,2,...,M). (3.39) 
fe=i 

The proof of (3.35) can be performed on the basis of the formulas (3.37)-(3.39). First, we 
multiply (3.37) by hji/pf and sum up with respect to i and j to obtain 

M ^ M M ^ 

^'M^j,M+iZp: = XI ^'i M+.i,k M+fiik -jo 
i,j=l i,j=lk=l 

MM ^ 

+ J] J]AM+,>M+i4'2i'2fc' {^,J = l,2,...,M). (3.40) 

ij=l fe=l ^» 



Note that for any function /, 



M 



E ^- = E E E (3.41) 

where the notation Xl^eSi implies that the dummy index /i runs over the set Si. Applying 
this rule to the left-hand side of (3.40) 

M ^ n ^ 

= E ^M+jM+i:^ ^ E E E ^M+i^M+t^^- (3.42) 

ij=l i,j=l ti&Si ueSj 

We modify L by taking into account the relations b^f^ — —bf^u and -D^^^ jy^^^ = D'm+ij.,m+i^ 
which follow from the skew-symmetry of the matrics D and D'. This leads to 



1 " / 1 1 \ 

^ 2 E E E ^'M+u,M+t, [ +-^)b„u 



1 



8 E E E ^M+.,M+M' (3-43) 

where in passing to the second line of (3.43), we used (3.14). It follows from (3.10b) and 
the formula (3.18) that 

X Xl^M+.-M+M= l^'(-ei;ej)l = |i?(-z,-ei;z,e,-)|, (3.44) 
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which, substituted in (3.43), gives 

1 " 

L = - ^ Cij\D{-z,-ei;z,ej)\. (3.45) 
On the other hand, using (3.38) and (3.39), the right-hand side of (3.40) reduces to 

M M 

We substitute the exphcit form of from (3.12) and take into account the symmetry 
= D'f^^, the first term of R is modified as 

^, / 1 1 \ Pi - Pa 



1 ^ / 1 2 1 \ 



It turns out by applying the formula (3.18) to (3.47) that 

M 

E ^^^^ = -|i^'(-z;z.O| - |i?'(-z.;zO| + -|i^'(-z..;z)| 

i,k=i Pi ^ 

1 1 
= 2'-^^"^' -z;z,z^^)| - |-D(-z, -z^;z,z^)| + 2l-^(~^' -Zrr; z,z) 

= -|L>(-z,-z^;z,zO|, (3.48) 

where in passing to the last line, we used the fact that any determinant which contains 
two identical rows (or columns) is zero. The similar procedure applied to the second term 
of R yields 

M 

J2D,,^ = \D{-z-Zrr)\. (3.49) 

i,k=l P"- 

Adding (3.48) and (3.49), we finally obtain 

R — |-D(— z; Ztt)\ — |-D(— z, —Zt'-i z, Zt-)|. (3.50) 

The desired relation (3.35) follows immediately from (3.40), (3.45) and (3.50), completing 
the proof of Eq. (3.7). 
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D. Remcirks 

1. Let C — {cij)i<ij<n be a real symmetric matrix whose dinagonal elements are zero, 
i.e., Cii — 0{i — 1, 2, n), and P — {Pij)i<i,j<n a regular matrix. Then, under appropriate 
orthogonal transformation Ui = Yl^j^iPiju'j, the quadratic form (1.7b) can be recast to a 
canonical form 

p q 

F = Y,uf-J2 ^'p+v {P + Q<n), (3.51) 

i=l i=l 

where p{q) is the number of positive (negative) eigenvalues of C, and p and q are deter- 
mined uniquely by C.^^ Note that since TrC = 0, p ^ and q ^ 0. Under the same 
transformation, the system of bilinear equations (3.6) and (3.7) can be converted into the 
system 

DyD,f-g[^fg[, (i = 1, 2, + g), (3.52) 



Dlf-f = 




(3.53) 



where u'^ = g[/ f{i = 1, 2, ■■■,p + q). For example, if Cij = 1 {i ^ j), ca = 0, then p = 1 and 
q = n — 1 since the eigenvalues of C are n — 1 (simple root) and —1 ((n — l)-plc root). 

2. When F is a positive definite quadratic form of {i — 1,2, ...,n), we can put p — n 
and g = in (3.52) and (3.53) provided that C has n distinct positive eigenvalues. The 
system corresponding to (1.7) becomes 

_ 1 

UijXt — 2 

If we consider the continuum limit n — >■ oo for (3.54), then we have a (2 + l)-dimensional 
nonlocal PDE of the form 




(3.55) 



This equation is an analog of the (2 + l)-dimensional nonlocal nonhnear Schrodinger 
equation 

POO 

iut — Uxx + 2u / u — u{x,z,t), (3.56) 




, {i = l,2,...,n). 



(3.54) 
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arising from a continuum limit of the multi-component nonlinear Schrodinger equation. ^^'^^ 
By means of the hodograph transformation 

dy — rdx + u^dz^ rdr, dt — dr, (3.57) 

we obtain the parametric representation of the solution for Eq. (3.55) 

u^^, (3.58) 

where / = f{y,T) and g — g{y, z,t) satisfy the system of bilinear equations 

1 f'^ 

DyDrf-g = fg, D'J ■ f = - J g'dz. (3.59) 

We will discuss the integrability of the equation (3.55) in a separate paper. 

3. The bilinear equation (3.7) takes the same form as that of a coupled modified Korcwcg- 
de Vries equations proposed in Ref. 18 where the proof of the multisoliton solution has 
been performed by lengthy calculations using various formulas of pfaffians. Here, we have 
provided a novel proof relying only on an elementary theory of determinants. 

4. The coupled PDEs proposed recently in Ref. 19 

Ui,xt ^Ui- ^ CjkiUj^^Uk - UjUk,3:)Ui, {l = 1, 2, u), (3.60) 
l<3<k<n 

where the coupling constants Cjk are skew-symmetric, are transformed to the following sys- 
tem of bilinear equations through the dependent variable transformations Ui — gi/f {i — 
1,2, ...,n) 

DM-gi^fgi, (i = l,2,...,n), (3.61) 
DxDtf ■ f = ^ CjkD^gj ■ gk. (3.62) 

l<j<k<n 

Recall that the bilinear equation (3.61) coincides with (3.6) if we replace the variables x 
and thj y and r, respectively. We conjecture that the multisoliton solution of Eqs. (3.61) 
and (3.62) will be given by (3.16) where the matrix Bm has the form 

b^v = -Cij — fie Si,ue Sj {ii,u = 1,2,...,M (/i ^ = l,2,...,n {ij^j)), 

(3.63) 
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in place of (3.14). Obviously, the corresponding tau-functions / and Qi satisfy Eq. (3.61) 
since its proof does not depend on the explicit form of Bm except that it is a skew- 
symmetric matrix with the constant elements. For the 2-component system, we have 
checked that Eqs. (3.61) and (3.62) exhibit the 2- and 3-soliton solutions, i.e., Mi = 
M2 = 2, Ml = M2 = 3. The proof of the general multisoliton solution will be reported 
elsewhere. 

IV. TWO-COMPONENT SYSTEM 

Here, we consider the two-component system (1.8) in detail. We first show the in- 
tegrability of the system by constructing a Lax pair and then present the multisoliton 
solution. We also discuss an integrable system (1.9) which is closely related to the system 
(1.8). 

A. Integrability 

For the system (1.8), the equations (3.2) and (3.3) corresponding to Eqs. (1.8) read 

where the first of these equations comes from the ^/-derivative of the second equation of 
(3.2) with F = uv. The system of equations (4.1) can be derived from the compatibility 
condition of the system of linear PDEs 

^fy = = (4.2a) 

with 

- = V^IQ (4.2, 

where A is a spectral parameter. Note in this expression that Xy — -y/l — UyVy. Indeed, it 
follows from the condition ^y-,- — ^-j-y that 

Ur-Vy + UV -VU = O, (4.3) 

which yields Eqs. (4.1). Using (2.1b), we can rewrite (4.2) in terms of the original 
variables x and t 

'^^ = m, ^t = V-^, (4.4a) 
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with 



.v^ -1 



~ _ 1 /O -u\ 

2 U j ^ 4A U -ij ^ 2 \uvu. 



\ f UV UVUx 

-uv 



(4.46) 



This is a Lax pair for the system of equations (1.8). Note that when u = v, (4.4) reduces 
to the Lax pair for the SP equation.^ One can apply the inverse scattering transform 
(1ST) method to estabhsh the complete integrability of the system (1.8). 



B. Multisoliton solution 
1. A^-soliton solution 

The parametric representation of the multisoliton solution of Eqs. (1.8) is given by 
(3.5) with n — 2 

ii=y, = (4.5) 

Here, we consider the case where both u and v contain solitons. Correspondingly, we 
set Ml = Ma = iV and M = 2N in (4.5). The tau-functions / and gi{i = 1,2) from (3.16) 
are represented by the following formulas: 



/ = Vf, 



^(^ = 1,2), 



(4.6a) 



with 



\D\ 



Gi = |i:>(-z, -ej;z,ej)| 





A2N 


I2N 








— I2N B2N 






^2N 


I2N 


0^ 




I2N 


B2N 0^ 


T 




-z 

























(i = l,2). 



Here, the 2N x 2N skew-symmetric matrices A2N and B2N have the elements 

A -(n \ n — ^' ~ n^i+ii — ^' ~Pi ^ ^ 
■^2N — [0'ij)l<i,j<2N, — ; e = ; ZjZj, 



Pi +Pj 



1 



= PiV + -r + ^io: (^ = l,2,...,2iV), 
Pi 



B. 



2N 



NxN B]\fxN 



~Bj^y^]\f Onxn 



(4.66) 



(4.6c) 



(4.6d) 
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Bnxn — {biN+j)i<i,j<N, biN+j — -7-^Y^^~^T^^ (^)i — 1) 2, A?"), (4.6e) 

4 Pi Pn+j 

and the 2A^-component vectors z and ei (i = 1, 2) are given respectively by 

^«^e^^...,e«-), ei = (l,l,...,l,0,0,^, 62 = (0, 0, 0, 1, 1, (4.6/) 



le^ - 

N N N N 



Note that the A^-soliton solution contains 4A'" complex- valued parameters pi, ^io (i = 
1,2,..., 2A'"). An alternative parametrization with the same number of the parameters is 
possible if one puts pAr+i = Pi {i — 1,2,...,N) and replace ^io and ^N+io by + Inoi 
and ^io + ln6i (i = 1,2, ...,N), respectively where and bi are new parameters. In the 
following, wc present a few examples of solutions and investigate their properties. 
2. 1-loop soliton solution 

We give two types of tau-functions described above: 

f T , 1 (PlP2)^ . . „ X 

f = ^ + 77 [ \2 ^1^2, 9i = zu 92 = Z2, (4.7a) 

4(Pi +P2Y 

/ = 1 + ^^^^^^ zl, gi^aizi, 5^2 = 6i2;i, (01,02, pi > 0). (4.76) 
The solution corresponding to (4.7b) is calculated from (3.5) to give 

= — 4 /p- sech(^i + 5i), = — ^/ — sech(^i + 5i), (4. 80) 

Pi V ^1 Pi V oi 

x = y--tanh(6 + 5i), 5i = In f ^^^^^^y (oi, 6i,pi > 0). (4.86) 



Pi 

A profile of m is depicted in Fig. 1. It represents a loop soliton with the amplitude ^\J^ 
and the velocity Ci = l/p\. Note that the amplitude of the loop soliton is defined by the 
maximum value of u which is attained at ,^1 = —5i in the present example. The property 
of V is the same as that of u except the amplitude given by ^-y/^- By comparing (2.16) 
and (4.8), we see that the loop soliton has the same structure as that of the loop soliton 
of the SP equation. 
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X 

FIG. 1. The profile of the 1-loop sohton solution u with the parameters pi — 1.0, oi = 0.5 
and bi — 1.0. 

3. 2-loop soliton solution 

As in the case of the 1-soliton solution, we write down two types of tau-functions for 
the 2-soliton solution: 




^ 1 {piP2P3PA)\pi-P2f{p?.~PAf z,z.,z,z^, (4.9a) 

16 (pi +P3)2(p2 +P3)^(Pl +P4)^(p2 +P4)^ ^ ^ ^ ^' 

- + + l{p.+v.?{v. + P.Y '''''' + 4 (^TT^IFb^T^ ^'-'^^ 

(4.10.) 

,1 , Pl{Pi-P2f 2,1 , pI{Pi-P2? 2 

5(1 = ai2;i + 02^2 + 7^010201^ \ ^1^2 + -rF^ia-^h—, ■ ZiZ^, (4.106) 

16 (Pi+P2)^ 16 {Pi+P2y 

h If. . ^ ^ u PljPi-P-i? 2 ,1 ^ ^ Pl{Pi-P2? 2 

5'2 = Oi^i + 022;2 + —aihih2—, ■ 2^2 + -77020162^ ^ ^1^2- (4.10c) 

16 {P1+P2Y 16 (Pi+P2)^ 
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We consider the 2-soliton solution corresponding to the tau-functions (4.10). Figure 2 
shows the time evolution of of the 2-soliton solution u. It represents the interaction of two 
loop solitons, each takes the form of the 1-loop soliton given by (4.8), as we demonstrate 
now. 




FIG. 2. The time evolution of the 2-loop sohton solution u with the parameters pi — 
0.5, p2 = 1.0, tti = 1.0, 02 = 2.0, bi = 1.0, 62 — 2.0 and Xio = X20 = 0. 

We investigate the asymptotic behavior of the solution u. To this end, we assume 
< pi < P2 and Qi > 0, bi > {i — 1,2). Then, an asymptotic analysis similar to that 
developed for the 2-loop sohton solution of the SP equation shows that as i — > —00, u 
behaves hke^ 

u = U1 + U2 ^ — . sech(^i + d[) + — . / ^sech(^2 + (^2), (4.11a) 
Pi \ bi p2 V ^2 

£1 2 2 4 

X + cit — xio tanh(^i + S[) , for ui, (4.116) 

Pi Pi Pi P2 

(^1 2 2 
X + C2t - X20 ^ tanh({2 + ^2) , for M2, (4.11c) 

P2 P2 P2 



where 



P 



,2 ' 



In 



Pi , = In 



y/oibi 



-Pi 



P1-P2 

Pi +P2 



, (z = l,2). {A.lld) 



As t ^ -|-oo, on the other hand 

2 [al 



u^Ui + U2^ — . I sech(^i + 5i) + — W ^ sech(^2 + So), (4.12a) 
Pi \ bi p2 V ^2 
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t 2 2 
X + cit — xio ~ tanh(^i + 5i) , for ui, (4.126) 

Pi Pi Pi 

2 2 4 
X + C2t — X20 tanh(^2 + ^2) , for ^2- (4.12c) 

P2 P2 P2 Pi 

We observe that the solution u sphts into two loop solitons as time evolves, each of which 
has the form of single loop soliton. The only effect due to the interaction of two loop 
sohtons is the phase shift. To see this, let Xic be the center position of the ith sohton. 
Then, it follows from the asymptotic forms (4.11) and (4.12) that 

xic + Cit - XiQ , X2c + C2t-X2o^ , [t ^ -oo), (4.13a) 

Pi Pi P2 P2 P2 

5i 2 5' 2 4 , , 
xic + cii - xio ~ , X2c + C2t - X20 . (t ^ +00). (4.136) 

Pi Pi P2 P2 Pi 

Since two solitons propagate to the left, the phase shift of the ith soliton can be defined 
by the relation 

Aj = Xi^{t -00) - Xic{t +00), (i = 1, 2). (4.14) 
Thus, from (4.13) and (4.14) one has 

A. = -ll„(!^l^)^l. (4.15a) 

Pi \Pl+P2j P2 

A. = llnf^^y + 1. (4.156) 

P2 VP1+P2/ Pi 

The same calculation can be applied to v as well. The corresponding asymptotic formulas 
are obtained simply by interchanging and 6^ {i — 1,2) in the above expressions. It 
should be remarked that the above formulas for the phase shift do not depend on a^ 
and 6j (i = 1, 2) and are determined only by the amplitude parameters pi and p2. They 
coincide with those of the 2-loop soliton solution of the SP equation.^ A novel feature of 
the solution in the present 2-component system is that the large soliton propagates slower 
than the small soliton if the inequality — < — holds. This fact is seen from the 

^ pi Y oi P2 Y 02 

asymptotic forms (4.11a) and (4.12a) of the solution with the velocities of Ui and U2 being 
given respectively by ci = 1/pf and C2 = l/p| (c2 < ci). 
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4. Breather solutions 

The breather solutions are constructed from the soliton solutions by following the same 
manipulation as that used for the soliton solutions of the SP equation.^ Here, we present 
the 1-breather solution. In this case, we put 

= a + i6 = p*, = A + i// = ^2*o> «i = CKie'*^' = ^'i = /^le'^' = b*, (4.16) 

in (4.10) to obtain the tau-functions /, gi and g2- Here, a,b,ai and f3i are positive 
constants. A, /i, 0i and ipi are real constants and the asterisk denotes complex conjugate. 
After a few calculations, we find the following expressions: 

with 

/ = 6^ cosh^(^ + ^o) + cos^(x + Xo + (^') - (a^ + h^) sin^ 5, (4.176) 
^1 = sin(xo-5)sin(x+Xo+5') cosh(6'+6'o)- cos(xo-5) cos(x+Xo+5') sinh(6'+6'o), (4.17c) 
^2 = sin(xo+(^) sin(x+Xo+(^') cosh(6l+6'o) - cos(xo+(^) cos(x+Xo+(^') sinh(6'+6lo), (4.17d) 
where 



0? + ^" ) ' ^ V a' + 6V 

e'° = ^V«i/5i(«' + &'), tanxo = ^, 5 = ^(0i-^i), 5' = ^(^i + V^i)- (4.186) 

Substituting (4.17) into (3.5), we obtain the parametric representation of the 1-breather 
solution: 

u = : — , V = : — , (4.iya) 

+ 62 / + 62 / ^ ^ 

2ab b sinh 2(g + gp) + a sin 2(x + Xo + ^^0 4a 
^"^ a2 + 62 / a2 + 62- ^^-'^^^ 

Both u and include two different phases 9 and x- The former characterizes the envelope 

of the breather whereas the latter governs the internal oscillation. Figure 3 shows the time 

evolution of the 1-breather solution u. It represents a oscillating localized pulse moving 

to the left. Contrary to single loop soliton, the profile of the pulse is nonstationary in the 
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comoving coordinate system. An inspection shows that the solution (4.19) exhibits singu- 
larities as encountered in the case of the breather solution of the SP equation. Therefore, 
certain condition must be imposed on the parameters a, b and 5 to produce the regular 
breather. However, we do not pursue the detail here. 




FIG. 3. The time evolution of the 1-breather solution u with the parameters a — 0.1, b — 

0. 5. ai = 1.0, pi = 2.0, 01 = 0, V^i = 7r/2, and A = // = 0. 

C. Related integrable system 

1. 2-component system 

The 2-component system (1.8) can be transformed to another integrable system (1.9) 
by a simple transformation. To show this, we put 

u = u + iv, v = u — w, (4.20) 

and substitute this into Eqs. (1.8), we obtain a system of equations for u and v 

Uxt = u + ^[{u^ + v^)u^]^, (4.21a) 

Vxt^v + ^[{u' + v'')vxl. (4.216) 

This system is a special case of the n-component system (3.54) with n = 2. 

2. A^-soliton solution 

The parametric representation of the A^-soliton solution for Eqs. (4.21) can be ex- 
pressed in the form 

u-f v^f - = y-24, (4.22) 
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where the tau- functions gi,g2 and / satisfy the system of bihnear equations 

DyDj-g^^rgi, (z = l,2), (4.23a) 

Dlf-f = ligi + gl)- (4.236) 

Here, we consider real- valued solutions u and v for the system (4.21). The tau-functions 
representing the A^-soliton solution are obtained from (4.6) by putting at = P*, Ci+No — 
{i — 1,2,...,N). Then, the expressions corresponding to (4.6d), (4.6e) and (4.6f) 
become 



^2 ^ij V Pi+Pj Jl<i,j<N \ Pi+Pj 



l<i,j<N 



Ki. _ 

(4.24a) 

z= (e^^e«^...,e«'^,e«^e^2,...,e«^). (4.24c) 

We can see that the tau-functions f,gi and g2 (4.6) with (4.24) satisfy the conditions 
f — f* and g2 — g*, which, combined with (4.5) and (4.20), gives 

/ = ^1 = ^(^1 + ^1*), ^2 = ^(^1-^!). (4.25) 

As in the case of the A^-soliton solution of the 2-component system (1.8) (see Sec. B), 
we have an alternative parametrization of the A'"-soliton solution. Namely, we replace 
by ^jo + ln(aj -|- ibj) {j — 1, 2, N) where aj and bj are real parameters, and then take 
Pj and {j — 1,2, ...,N) be real in the expressions of the tau-functions. This procedure 
yields the tau-functions corresponding to (4.7b) and (4.10), for example. Actually, for the 
1-soliton solution, the corresponding tau-functions are given by 

f ^i + ^ial + bDplzl, gi^aizi, g2^biZi, (4.26) 
and for the 2-soliton solution, they read 

/ = 1 + ^(«? + bDpjzf + ^{a,a2 + bib2)j^^^^^ ^1^2 + ^(4 + bDpjzj 
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256 ^""^ ^''^''2 (p, + p,)4 ^^1^2) , (4.27a) 

5(1 = ai2;i + 02-2:2 + — 02(01 + 61)— 2:2 + — 01(02 + 62)— — Ziz^, (4.276) 



16 (pi+P2)^ 16 (P1+P2) 



> , 1 , / 2 , ,2nP? Pi -P2 2 , 1 , / 2 , ,2NPi Pi -P2 2 07 ^ 

^2 = 61^1 + 62^2 + 7^62(01 + 61)— ■ 7^^1^2 + 7^61(02 + 62)— ■ ^ziz^. (4.27c) 

16 (Pl+P2)^ 16 {Pl+P2)^ 

It can be seen that substitution of (4.26) and (4.27) into (4.22) produces the 1- and 2-loop 
sohton solutions, respectively. 
3. Breather solutions 

One can confirm by direct calculation that the tau-functions (4.27) satisfy the bilinear 
equations (4.23). In the process, the reality of the parameters has not been used. This fact 
enables us to extend the range of the parameters to complex values. Thus, the breather 
solutions are constructed from the soliton solutions by applying the procedure developed 
in Sec. IV B. In practice, according to the parametrization (4.16), one has 

/ = 4 ^2(.-^.o)/, 'JT, 9u ~92- 92, (4.28a) 

with 

/ = 6' cosh'(e + ^1) + a' cos'(x + Xo + «) + ^(o' + 6^) (^^^^ - , (4.286) 

gi = sin(xo-0i+«) sin(x+Xo+'«) cosh(^+^i)- cos(xo-0i+'«) cos(x+Xo+'«) sinh(^+^i), 

(4.28c) 

§2 = sin(xo-V'i+«) sin(x+Xo+«) cosh(6l+6li)- cos(xo-V'i+«) cos(x+Xo+«) sinh(6'+6li), 

(4.28d) 

where the parameters 9i , 7 and k are defined by 

e^i = A v'^2Tp^ ^ _ [q,4 ^ 2alPl cos 2(01 (4.28e) 

al sin 201 + /Sf sin 2z/'i 
tan 2k = — -75 -— , (4.28/) 

al cos 201 + (51 cos 20i ' \ J ' 

and d,x ^"^^ ^0 are already given respectively by (4.18a) and (4.18b). Substituting (4.28) 
into (4.22), we obtain the parametric representation of the 1-breather solution: 

i^^^a. « = ^4|=f. (4.29a) 

7Va'' + 6^ / 7Vo'' + 6^ / 

28 



x = y- 



2ab b sinh 2{0 + 0i) + a sin 2(x + Xo + 4a 



0,2 + 52 



/ 



(4.296) 



Of particular interest is a circularly polarized wave for which the solution exhibits a 
simple structure, as we shall now demonstrate. In this case, we put ai — f3i, Xo~<Pi+i<' — f 
and 01 — ■^i = I to obtain the tau-functions 



4a2 



e^^ ^1 = 2Q;ie'cos(x + 0i), ^2 = 2q;i e' sin(x + 0i). (4.30) 



Then, the solution takes the form 



u — 



2a cos(x + 01; 
a2 + 62 cosh(^ + 9'o 



a2 + 62 cosh(^ + Q'q 



2a 



2a 



X = y 



a^ + b"" 



tanh(^ + O'q) - 



2a 



+ 62' 



(4.316) 



The parametric solution (4.31) represents a nonsingular breather if the inequality < 
a/6 < 1 holds. Figure 4 shows the time evolution of u{= u) given by (4.31). 




FIG. 4. The time evolution of the 1-breather solution u{= u) with the parameters a — 
0.1, 6 = 0.5, ai ^ 1.0, 01 = and A = 0. 

One can show by a direct calculation that the solution (4.31) satisfies the integral 
relations 

/oo poo 
udx^O, / vdx^O, (4.32) 
-00 J —00 

implying that both u and v are zero mean fields. This fact indicates clearly an oscillating 
character of the solution. Note that the above relations represent the conservation laws 



29 



derived from the system of equations (4.21) for localized waves. It is interesting to observe 
that in the small amplitude limit a/h — > 0, the profile of u bears resemblance to that of 
the soliton solution of the nonlinear Schrodinger equation. 

D. RemEirks 

1. The system of equations (4.1) is equivalent to a coupled dispersionless system for the 
variables r = r(x, t), s — s{x, t) and q — q{x, t) 

Qxt + {rs)a: = 0, r^t - '^Qxr = 0, s^t - '^QxS = 0. (4.33) 

The Lax pair associated with the system (4.33) has been obtained and the 1ST has been 
applied to it to construct soliton solutions. In particular, 1- and 2-soliton solutions have 
been presented for r, s and q^- Here, we present the formula for the general multisoliton 
solution for the first time. 

2. The system of bilinear equations (4.23) can be derived from the system (4.33) with a 
reduction s = r* through appropriate dependent variable transformations. ^-"^ See also an 
analysis by means of the IST.^^ 

V. CONCLUSION 

In this paper, we proposed a novel multi-component system associated with the SP 
equation and constructed its multisoliton solutions in terms of pfaffians. We also consid- 
ered the equations reduced from our system. In particular, the 2-component system (1.8) 
was found to be completely integrable for which the explicit Lax pair was presented. We 
also provided the loop soliton and breather solutions for the system and investigated their 
properties. We also addressed the system (1.9) which stems from the system (1.8) by a 
simple transformation. In conclusion, we shall discuss some open problems associated 
with the multi-component system under consideration. 

1. One interesting issue to be resolved in a future work is the proof of the complete 
integrability of the n-component system (1.7) by using the 1ST. To construct the Lax pair 
for the system, one way will be to start from the system of bilinear equations (3.6) and 
(3.7) to obtain the Backlund transformation among the tau- functions and then derive the 
scheme of the 1ST following the standard procedure in the bilinear formalism. 
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2. Other issues to be reserved for detailed study have aheady been described in Sec. 
HID. Of particular importance is the construction of the multisoliton solution of the 
n-component system (3.54) with n > 3. Unlike the 2-component system, the linear 
transformation such as (4.20) does not exist to convert the system (1.7) to the system 
(3.54). Hence, one must solve the system of equations (3.52) and (3.53) with p = n and 
q = 0. It will be a relatively simple task to obtain the 1- and 2-soliton solutions analogous 
to (4.26) and (4.27). Nevertheless, a systematic approach is necessary to construct general 
multisoliton solutions. 

3. The system of equations (1.5) has been derived as a unidirectional model describ- 
ing the propagation of circularly polarized ultra-short pulses in a Kerr medium. The 
solution of breather type has been obtained by means of an analysis as well as numer- 
ical computations. ^^'^^ However, in view of the extremely complicated structure of the 
breather solution, it seems to be unlikely that the system admits multibreather solu- 
tions as well. Thus, we suspect the complete integr ability of the system even if it has 
passed the Painleve test. On the other hand, although the difference between (1.5) and 
(1.9) is the location of the x deivative on the right-hand side, the latter shares many 
common features to the integrable systems such as the complete integrability and the 
exsistence of multisoliton solutions. At present, however, the relevance of the system to 
the description of the dynamics of ultra-short pulses in optical fibers is not clear. Nev- 
ertheless, it would be of interest to examine the possibility of the system as a physical 
model for the two-component generalization of the SP equation. 
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